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I. INTRODUCTION 
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The nucleus-nucleus interaction potential has been studied within the relativistic mean field the- 
ory. The systematics of the relativistic effects have been investigated by analyzing the relation 
between the potential and the bombarding energy as a function of the impact parameter. It is 
shown that the potential barriers are noticeably sensitive to the bombarding energy for a given 
impact parameter. At large bombarding energies the slope at the potential edge decreases with 
the impact parameter. Comparisons with a non-relativistic treatment shows that relativistic effects 
cannot be ignored at bombarding energies around and larger than 100 MeV/nucleon. 
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A consistent theoretical treatment of relativistic many-body systems is still a challenge in many areas of physics. 
The effects of the Dirac sea imply corrections due to vacuum fluctuations and other features of quantum field theory 
may also become manifest. An important aspect of special relativity, namely causality, is only verified if sea effects are 
properly incorporated. Unfortunately, up to now it is proved to be very difficult to find unambiguous experimental 
evidence for such effects. That might not be the case for static properties of nuclear matter for which a successful 
project carried out during many years has shown that nuclear energies for the ground and lower excited states of 
nuclei are well described by a relativistic mean- field treatment of the collection of nucleons [H-@ . 

Almost nothing has been done to understand the rather relevant problem of the influence of relativity on many- 
body scattering. Usually, this is accounted quite straightforwardly in quantum field theory for two-body scattering 
including the sea effects with help of perturbation theory. But the scattering of relativistic composite objects where the 
compositeness becomes an active part of the scattering process is still a very difficult theoretical problem. This is an 

j. obviously important problem in nuclear physics, due to an inherent theoretical difficulty in defining a nuclear potential 

, between many-body relativistic systems. Nonetheless, nuclear potentials are often used in the analysis of nucleus- 
es ■ nucleus scattering at energies of 100 MeV/nucleon and higher. In fact, many rare isotope facilities use reactions at 
these energies in which relativistic effects are expected to be of the order of 10% or more for some processes of interest 
' ■ For proton- nucleus scattering a successful approach, known as "Dirac phenomenology" , has been use for many 

\ years [l(|. But for nucleus- nucleus collisions a reasonable account of these features has not yet been accomplished. 
t-H . Due to retardation, an attempt to use a microscopic description for nucleus-nucleus potentials starting from binary 
collisions of the constituents is not possible as a nucleus- nucleus potential requires a simultaneous interaction between 
all constituents. In the case of nucleus-nucleus collisions a similar approach as in the Dirac-phcnomcnology case can 
be followed up by using a mean field theory. 

The nucleus- nucleus potential used in the analysis of elastic scattering experiments is a very complex object which 
contains information about all possible inelastic processes. The presence of inelastic channels lead to, among other 
things, the introduction of an imaginary part for the nucleus-nucleus potential. This potential is usually termed 
"optical potential" . Although several microscopic models have been developed to calculate the optical potential, due 
to its complexity, it is often parametrized by a set of phenomenological functions depending on the distance between 
the center of mass of the interacting nuclei. The parameters set (often with ambiguities) is obtained by fitting elastic 
scattering data whenever available. Needless to say that for many systems (e.g. most reactions with radioactive beams) 
the elastic scattering data are unavailable and one has to resort to theoretical constructions of the optical potential. 
A popular method is to build a folding potential in which the real part of the optical potential is obtained from an 
integral over the nucleon-nucleon densities weighted by their individual (sometimes density dependent) interactions. 
The nuclear densities arc taken as frozen, with all nucleons interacting simultaneously. A method to deduce the 
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magnitude of the effects of relativity in nucleus-nucleus interactions presented in Rcfs. [7H2| where the nucleus- 
nucleus potential for a relativistic projectile was modified in a similar way as the scalar part of the electromagnetic 
interaction. This procedure has shown that relativistic corrections lead to appreciable changes of inelastic processes 
involving 100 MeV/nucleon projectiles. 

In this work, we follow the steps of relativistic mean field theories to construct a relativistic nucleus- nucleus potential. 
Our goal is to study how relativistic effects would be manifested in the coordinate dependence of a nucleus-nucleus 
potential. We use a relativistic mean field theory for each nucleus separately for consistency. The effect we are looking 
for is obtained by boosting the projectile motion to the target frame of reference (or to the laboratory frame of reference, 
if both projectile and target are boosted) to analyze how the interaction energy is modified. The interaction energy is 
defined below as the total energy of the system when the nuclei are at a given finite separation subtracting from it the 
sum of their individual energies when they are infinitely separated. This is an oversimplification of the problem and 
it also only allow us to obtain the real part of the interaction potential. In order to study the influence of relativistic 
effects on elastic cross sections we use, for simplicity, an imaginary part of the optical potential having the same radial 
dependence as the real part. 



II. ENERGY DENSITY FUNCTIONAL FOR SINGLE AND TWO-NUCLEUS SYSTEMS 

Based on the meson exchange theory, finite nuclear systems can be well described by an energy density functional 
within the relativistic mean field approach [H, 0. In the mean field approach, the nucleons are treated as point- 
like particles interacting by the exchange of mesons and photons. The energy functional [IMsl for a single nucleus, 
associated with a-, u>-, p-mesons and photon (A) exchange, is given by (h = c = 1) 

E = Jdr^a (-H ■ V + M) tp a + ^ J2 J drdr ' E Mr)Mr')Fo(r, r')D+(r - r')A(r')i> a (r), (1) 

a <p—(7,uj,p,A ab 

where the two-body interacting matrices read as 

r<r(r,r')=-g (r (r)g (7 (r'), (2a) 
r u (r,r') = + (ft/Or (ft-7/*) P / . (2b) 
r„(r, r') = + {g P i^f) r ■ {g P i P f) r , , (2c) 

TA(r, r')=+ e - [f (1 - r z )] r ( 7 „ (1 - r 2 )) r , . (2d) 

In the energy functional (fT]), tp a denotes stationary single particle states, M is the rest mass of nucleon, and g a , g u , 
and g p are the coupling constants with baryonic density dependence pH - |l5j . The propagators of mesons (<f) — a, uj, 
and p and for the meson mass) and photon (A) are given by 

1 e ™A r - r '\ 1 

D 1> = -a i 7T> Da= \ n- ^ 

in \r — r'| \r — r'\ 

In this paper, we use bold types to denote the space vectors and arrows for isospin vectors. Based on the energy 
functional (fT]), the single particle configurations (ip a ) can be obtained from a self-consistent iterative procedure, 
aiming at describing the physical properties for the ground state of the nucleus such as binding energy, radii, density 
distributions, etc. 0-0]. 

The relativistic energy functional (fT]) described above for a single nucleus can be also extended to describe the 
interaction between two nuclei, i.e., the interaction between a target and a boosted projectile. Here we only consider 
same target and projectile nuclei. As shown in Fig. fT] the two-nucleus system consists of a target and a projectile with 
velocity v, separated by a distance R. For simplicity, we assume straight-line dynamics, appropriated for high-energy 
collisions. A given point in space is denoted by Tt = (xt,yt, Zt) in the frame of reference of the target with respect to 
its center of mass. For a collision with impact parameter b this point can be related to its coordinate r p = (x p , y p , z p ) 
in the frame of reference of the projectile by means of a Lienard-Wiechert transformation. If v is taken as the z axis, 
this transformation reads 

x p = x t + b, x t = x p — b, 

Up = Vt, Vt = y P , 

z p = "f{z t + RcosO), zt= -f(z p — RcosO), (4) 



3 




FIG. 1: (Color online) Two-nucleus system containing rest target and boosted projectile with the velocity v. The ellipse 
represents the boosted projectile density as seen from the target reference frame. R and b denote the distance and impact 
parameter between two nuclei in the target reference frame. 



where the impact parameter is related to the relative distance between the center of mass of the nuclei by b = i?sin#. 
The Lorentz factor is 7 = (1 — and /3 = v/c. 

The total energy functional for two- nucleus system shown in Fig. |T| is given by 

E(A t ,A p ,v) = E(A t )+E{A p ,v)+E(A t ,A p ,v), (5) 

where E(A t ) and E(A p ,v) denote respectively the target (with mass number A t ) and the projectile (with mass 
number A p ) energy functional, and E(A tl A p ,v) represents the interacting potential energy. For the isolated target 
and projectile, the energy functional E(A t ) and E(A p ,v) = E(A p ) can be determined from Eq. ((T|). Similar as the 
second term in right hand of Eq. JT]), the interacting one E(A t ,A p , v) can be obtained as 

E(A t ,A p ,v) = J ' dr j ^^'E^( r )^ b ( r ') r ^ r ' r ')^( r - r ')^ fc ( r ')^( r )' ( 6 ) 

4)—a,u,p,A ab 

where the integral over r and r' are done within the same reference frame, e.g., the target inertial frame. The single 
particle configurations ipt a and ip p & (subindex t for target and p for projectile) can be determined from the individual 
self-consistent calculations in the target and projectile frames, respectively. Because the collisions arc fast, we adopt 
a frozen configuration scheme in which the rearrangement effects from the nucleus-nucleus interaction are neglected 
in determining the single configurations of each nucleus. 

Taking a and uj as example, the corresponding contributions in Eq. ([5]) can be written as, 

E„=-- dr t dr' p g a {r t )p Stt (r t )D a (r t -r' t )p s , p {r' p )g a (r' p ), (7) 
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E u = + J dr t I ' dr' p g u (r t )p btt (r t )D u (r t - r' t )p biP {r' p )g u {r' p ). (8) 
The scalar and baryonic densities p s and pb can be obtained as 

Ps(r) =^2$ a {r)4> a (r), Pb(r) = ^^^hVaW, (9) 



where the sums are restricted within one nucleus. 

In the interacting energy functional ©, there exist two types of couplings, the scalar (cr) and the vectors (w, 
p, and ^4). For a straight-line motion, the Lorentz boots are: {pb, P (f p ), 0} — > { r ypb, P {r p ),'ypb :P (r p )(3}. The scalar 
coupling remains invariant under the Lorentz transformations because dr' t = dr' p /j and the Lorentz contraction of 
the longitudinal distance cancels out the Lorentz enhancement of the field. This contrasts with the vector coupling 
channels where the Lorentz transformation leads to a strong dependence on the projectile velocity. 
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In energy functional (j6)), the coupling constants are functions of density p v = j^j^ and = ipt^yPt + "^pl^^p, 
and the Lorentz transformation from the target to projectile, and the reverse, leads to 

Pv{r t ) =\JPb,t( r t) + ^lPb,t{r t )pb, P {r p ) + p£ iP {r p ), (10) 

Pv(r p ) =^plt(rt) + 2jPb,t{ r t)PbA r p) + Pl P ( r p)- (H) 

The method described above allows us to calculate the interaction potential of the two nuclei in terms of the single- 
particle states calculated within a relativistic mean field theory. In practice, the calculations are very complicated 
because the advantage of treating wavefunctions and densities in terms of angular momentum expansions in spherical 
basis is lost. Due to Eq. (TJ|, the boosted quantities also include Lorentz 7 factors within the polar angles of the 
single-particle states for the projectile rendering a very complicated description of their radial and angular dependence. 



III. RESULTS AND DISCUSSIONS 



In this work, we take 12 C- 12 C as a typical example of a two-nucleus system in Fig. [T]and study the nucleus-nucleus 
interacting potential energy according to the method described above. In the first step, each nucleus is treated in their 
reference frame independently to obtain the single particle configurations, from which the nucleus-nucleus interacting 
potential energy is determined via Eq. ©. The effective interaction PKDD is utilized to determine the single 
particle configurations as well as the nucleus-nucleus potential. 

In Fig. [2] we show the nucleus-nucleus potential energies V(R) = E [sec Eq. ©] as a function of the distance 
R for different values of Lorentz factor 7. The impact parameter is set as b = fm, and only the surface part of 
the potential has any useful application. The inset gives the results from R = 5 to 10 fm with enlarged scale. The 
nucleus-nucleus potential gets contributions from different fields at different times. First the nuclei approach from 
a fairly large distance, there exist very little overlap between the target and projectile densities and the potential is 
mainly contributed by the long range interaction - the photon field (A), which yields the Coulomb barrier. As the 
distance decreases to, e.g., R ~ 7.5 fm for 7 = 1.0, the densities start to overlap and the nucleus-nucleus potential 
bend toward a negative value. This is due the contribution of the meson fields leading to attraction. When R < 5 fm, 
the nucleus- nucleus potential have only net attractions. For the chosen nucleus 12 C, the self-consistent calculation 
with PKDD gives the neutron, proton and total radii, respectively, as r n = 2.26 fm, r p = 2.28 fm, and r — 2.27 
fm. These values (multiplied by 2) are consistent with the evolution of the nucleus-nucleus potential with respect to 
separation distance R, as seen in figure Fig. [5J 




R (fm) 

FIG. 2: (Color online) The interaction potentials (MeV) as functions of the distance R (fm) between rest target and boosted 
projectile nuclei for impact parameter b = fm and different values of 7. The inset presents the potentials from R =5 fm to 10 
fm with an enlarged scale. 

The nucleus-nucleus potentials shown in Fig. [2] also display a systematical behavior with respect to the Lorentz 
factor 7 which represents the boosted energy of the projectile, roughly (7 — l)M per nuclcon. With increasing 
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bombarding energy, the nucleus- nucleus potential well becomes shallower, as seen in Fig. [2j The point where the 
potential becomes attractive does not change appreciably with the bombarding energy. But, as shown in the inset, 
the potential barrier does have a strong dependence on the bombarding energy. This hints to a destructive additional 
(because of the boost) cancellation between the attractive meson fields and the repulsive photon field. These additional 
cancellations also tend to lead to a shallower potential. The filled circles in the inset show the tendency of this effect 
to stretch the nucleus- nucleus potential outwards as the energy increases. 

In order to have better understanding about the systematics of the nucleus-nucleus potentials, we show in Fig. [3Jthe 
detailed contributions from four coupling channels, (a) w-vector, (b) cr-scalar, (c) Coulomb vector, and (d) p-vector 
couplings. The corresponding insets present the results from R = 5 to 10 fm with enlarged scales. In Fig. [3] (d), 
it is clearly seen that the contribution from the p field is very tiny since in 12 C the nucleons are isospin saturated. 
Among the contributions to the nucleus-nucleus potential, the attraction mainly originate from the a scalar coupling, 
whereas the repulsion is mainly due to the uj and Coulomb fields. Starting from a fairly large separation distance, the 
sole contribution photon field yields and increasing Coulomb repulsion. At R ~ 10 fm the a and u> fields start to act 
with distinct contributions to turn the Coulomb repulsion into nuclear attraction. From Fig. 02a) and (b), one sees 
that the w and a contributions increase rapidly when the target and projectile approach. This is mainly due to the 
short range character of the meson exchange interaction. 

I — ' — I — 1 — I — 1 — I — ' — I — 1 — I — ' — I — ' — I — ' — I — ' — I I — ' — I — 1 — I — ' — I — ' — I — ' — I — ' — I — ' — I — ' — I — 1 — I 




R (fm) R (fm) 

FIG. 3: (color online) Contributions to the interaction potentials from the u (a), a (b), Coulomb (c), and p (d) fields are shown 
as functions of distance R (fm) with different values of 7 for the same two-nucleus system shown in Fig. [I] The insets show 
the details from R = 5 fm to 10 fm with enlarged scales. 

As seen from the inset of Fig. [3] (c), the photon field has a rather weak dependence on the boosted energy whereas 
a stronger dependence is found from the contributions from a and w fields. As shown in the insets of Fig. |3j (a) and 
(b), the repulsive (w field) and attractive (er field) increase with the boosted energy considerably When the target 
and projectile approach a close distance, the curves with different bombarding energies tend to cross each another at 
some point, denoted by red circles in Fig. [3] In fact these are turning points, where the trend of the bombarding 
energy dependence is just reversed. Similar turning points can be also found in Fig. [2j where the nucleus-nucleus 
potential becomes attractive. 

In Fig. [5] and Fig. [31 the impact parameter b is set to zero. In Fig. @]we show the nucleus-nucleus potential as a 
function of the separation distance R with the impact parameter b = 0, 1, 2,3, and 4 fm. For comparison, the results 
with four bombarding energies are shown in Fig. |4l (a) 7 = 0, (b) 7 = 1.1, (c) 7 = 1.3, and (d) 7 = 1.6. As seen 
from Fig. E](a) an d (b), the nucleus-nucleus potentials with lower bombarding energy are not sensitive to the impact 
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parameter. With increasing bombarding energy, the potentials show a strong dependence on the impact parameter. 
As shown in Fig. 0] (c) and (d), the slope of the nucleus- nucleus potential at large distances becomes smaller with 
increasing impact parameter. It is also demonstrated by the insets in Fig. |4] 
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FIG. 4: (Color online) The nucleus-nucleus interaction potentials (in MeV) are shown as function of the distance R between 
the target and projectile nuclei for different values of 7 [(a) 7 = 1.0; (b) 7 = 1.1; (c) 7 = 1.3; and (d) 7 = 1.6] and different 
impact paramters 6 (fm). The insets show the results from R — 3.5 fm to 8 fm with enlarged scales. 



In order to asses the rclativistic modifications of the nucleus-nucleus potential on scattering observables, we consider 
the elastic scattering cross section of spherically symmetric nuclei. At high energies the scattering amplitude is well 
described by the eikonal approximation, 



/0?) = fc(4)+ik / Jo(qb)exp[ixc(b)]{l-cxp[i X N(b)]}bdb 



(12) 



where Jo is the Bessel function of zero-th order, q = 2/csin(t9/2) is the momentum transfer for the scattering angle -d 
and projectile momentum hk = fiv. The reduced mass is /1 = 77710/(1 + 7), where mo is the rest mass of 12 C, and 
v 



C\J (7 — l)/7- The Coulomb-scattering amplitude reads 

Z x Z 2 e 2 



2fiv 2 sin 2 (tf/2) 



exp {-w?ln [sin 2 (tf/2)] + in + 2i(f> } , 



(13) 



where r\ = 2Z\Z2e 2 /Hv, and <fio = argT (1 + r)/2). 

The eikonal phase xn and xc in Eq. (fT2|) are given by 



1 [°° 

XJv(fr)=--^y Uopt(b, z )dz , 



xc(b)= 2 -^^Hkb), 

hv 



(14) 



where xc{b) represents the Coulomb eikonal phase. 

As mentioned before, the imaginary part of the optical potential has to be included by hand. In terms of the 
real part of the nucleus-nucleus potential V(r) (r = (b, z) and r = \Jb 2 + z 2 ) calculated above, the nuclear optical 
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potential will be written as 



U opt (r) = V(r)+i\V(r). 



(15) 



The coefficient A is set arbitrarily to 0.8 in the calculations. The elastic cross section is then given by da/dfl = |/(t9)| 2 . 
As usual, the elastic cross sections will be described in units of the the Rutherford differential cross section 



% = |/Ruth.W| 2 = 



1 



Z^e 2 
2fj,v 2 J sin 4 (tf/2)' 



(16) 



Our results for 12 C+ 12 C at 7 = 1.05 (approximately 50 MeV/nucleon) and 7 = 1.3 (approximately 300 
MeV/nucleon) are shown in Fig. [5] (solid lines). We compare the cross sections with those calculated with non- 
relativistic dynamics are shown in Fig. [5] by dashed lines. On observes that the minima of the elastic cross sections 
are shifted to smaller values of the scattering angle. This is due to the smaller effective radius of the nucleus-nucleus 
potential when rclativistic effects arc included, as observed in Fig. [H-Fig. F or a discussion of the effects of size and 
diffuscness of optical potentials on scattering observables, see, e.g., Ref. |l6|. The diffusencss of the potentials also 
increase with the inclusion of rclativistic effects. This is shown manifest in Fig. Fig. [5] by a (slightly) faster decrease 
of the cross section as a function of angle. 
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FIG. 5: (Color online) Elastic cross sections da/dO, = |/($)| 2 in units of Rutherford cross sections |/Ruth. (#)| as functions of 
the scattering angle with different values of Lorentz factor 7. The dashed lines represent the results with non-relativistic 
reduction. 



In summary, we have investigated the effects of relativity on nucleus-nucleus potentials within the rclativistic mean 
field theory. The rclativistic effects have been studied by analyzing the dependence of the potentials upon bombarding 
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energies and impact parameters. It is found that for a given impact parameter the Coulomb barrier is softened with 
increasing bombarding energy. For large bombarding energies, the potential edge becomes increasingly flat with 
increasing impact parameters which indicates that the target and projectile have to get closer to compensate for more 
attraction. The studies carried out here are exploratory in the sense that a consistent theory of a nucleus-nucleus 
potential including relativity requires a much more elaborated theory. Such a theory is missing in the literature. Many 
of nucleus-nucleus experiments at energies of 100 McV/nucleon and above are presently been carried out around the 
world with the goal to extract spectroscopic information on rare nuclear species [ill ]. Our work shows that some 
modifications on the values of the extracted spectroscopic quantities might occur due to the relativistic dynamics 
missing in most methods used to construct a nucleus-nucleus potential, such as the folding models. 
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